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We study the directional properties of a radiation field emitted by a geometrically small system composed 
of two identical two-level emitters located at short distances and driven by an optical vortex beam, a Laguerre- 
Gaussian beam which possesses a structured phase and amplitude. We find that the system may operate as 
a nanoantenna for controlled and tunable directional emission. Polar diagrams of the radiation intensity are 
presented showing that a constant phase or amplitude difference at the positions of the emitters plays an essential 
role in the directivity of the emission. We find that the radiation patterns may differ dramatically for different 
phase and amplitude differences at the positions of the emitters. As a reults the system may operate as a two- or 
one-sided nanoantenna. In particular, a two-sided highly focused directional emission can be achieved when the 
emitters experience the same amplitude and a constant phase difference of the driving field. We find a general 
directional property of the emitted field that when the phase differences at the positions of the emitters equal 
an even multiple of 7r/4, the system behaves as a two-sided antenna. When the phase difference equals an odd 
multiple of vr/d, the system behaves as an one-sided antenna. The case when the emitters experience the same 
phase but different amplitudes of the driving field is also considered and it is found that the effect of different 
amplitudes is to cause the system to behave as a uni-directional antenna radiating along the interatomic axis. 

PACS numbers: 42.25.Hz, 42.25.Kb, 42.50.Gy 


I. INTRODUCTION 

The control of emission by nanoscale-size objects such as 
single atoms, quantum dots, dimers and oriented semicon¬ 
ductor polymer nanostructures is fundamentally important for 
applications in nanoscale optical manipulation, optical sens¬ 
ing, information processing and quantum communication [1- 
3]. The recent advances in nanophotonics have stimulated 
a series of experimental and theoretical works demonstrating 
that composite nano-systems can serve as nanoantennas [4-7]. 
Examples include structures composed of nanodimers, high- 
permittivity dielectric particles, metal particles, and atomic 
chains [8-13]. It has been demonstrated that these nano¬ 
systems may squeeze light into nanoscale volumes [14], en¬ 
hance the excitation and emission rate of individual emit¬ 
ters [15-17], and tune the luminescence spectrum [18], and 
the polarization [19]. Particularly interesting is the ability to 
control and tune the radiation pattern of a nanoantenna. 

A related interesting problem is the ability of a composite 
structure of nanoparticles to work as a frequency filter which 
could route different frequencies of an incident beam into dif¬ 
ferent directions. This feature has been demonstrated exper¬ 
imentally by Shegai et al. [20]. In the experiment, a pair of 
metallic nanoparticles, gold and silver, was deposited on glass 
at a very small distance. When illuminated with white light, 
the system scattered the red and blue components of the inci¬ 
dent light into opposite directions. 

Recently, we have developed a theory of directional emis¬ 
sion for the somewhat related problem of directional light 
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scattering by a system composed of two two-level atoms [21]. 
We have shown that the system can operate as a directional 
nanonantenna provided that the atoms are not identical with 
unequal transition dipole moments or different transition fre¬ 
quencies. We have found that a difference between the transi¬ 
tion dipole moments or between the transition frequencies of 
the atoms creates a phase shift between the dipole moments 
of the atoms which then leads to the directional light scat¬ 
tering. Thus, the crucial factor for the directionality is to use 
emitters of unequal dipole moments or different transition fre¬ 
quencies. In the experiment of Shegai et al. [20] the required 
phase difference between the dipole moments was achieved 
by using two nonidentical metallic nanoparticles of different 
plasmonic frequencies. It has also been shown that a metal- 
dielectric structure composed of a pair of dielectric and metal 
nanoparticles can lead to a directional light scattering [22]. 
These studies show that a constant phase shift between dipole 
moments can be achieved with nonidentical nanoparticles. 

However, a constant phase shift between two oscillating 
dipoles could be achieved with identical nanoparticles. It is 
the purpose of this paper to demonstrate that a controlled di¬ 
rectional emission can also be achieved in a system composed 
of two identical emitters. As we shall see below, it requires 
a driving field which could create a constant phase difference 
between the atomic dipoles or a constant difference between 
the field amplitudes at the positions of the atoms. We con¬ 
sider a system composed of two identical emitters located at 
short distances and driven by a laser beam with a structured 
phase and amplitude, an optical vortex beam like a Laguerre- 
Gaussian (LG) or Bessel beam [23]. The use of an LG or 
Bessel beam can ensure either equal amplitudes and a constant 
phase difference or equal phases and a constant amplitude dif¬ 
ference at the positions of the atoms. We are particularly in- 
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terested in the directionality of the emission that is produced 
by the LG beam which is applied in the following manner; (i) 
a constant phase difference is produced at the positions of the 
atoms, and (ii) a constant amplitude difference is produced at 
the position of the atoms. We demonstrate that a controlled di¬ 
rectionality of the emission can be achieved under these spe¬ 
cific driving configurations. Especially, the atoms placed at 
short distances may operate as a highly directional nanoan¬ 
tenna. Depending on the phase or amplitude differences at the 
locations of the emitters, the system may operate as a two- or 
one-sided nanoantenna. A simple physical interpretation of 
the sources of the two- and one-sided emissions is given in 
terms of the collective states of the system. 

The paper is organized as follows. In Sec. II, we review the 
basic properties of the laser field of a structured complex am¬ 
plitude, an LG beam, and provide a simple explanation of how 
the beam could be applied to create a constant phase and/or 
amplitude difference between two dipoles located at different 
points. In Sec. Ill we describe the model in detail and discuss 
the method we use to to determine the radiation pattern of the 
system. The general expression for the radiation intensity is 
presented along with a brief discussion of its directional prop¬ 
erties. In Sec. IV we give illustrative figures of the direction¬ 
ality of the emission for two conhgurations of the LG beam 
at the positions of the atoms, either the same amplitudes and 
a constant phase difference or the same phases and a differ¬ 
ence between the amplitudes. We discuss different cases and 
present a detailed analysis of the directional properties of the 
radiation pattern. Finally, in Sec. V, we summarize our results 
and conclude. 


II. AMPLITUDE AND PHASE DIEEERENCE IN A 
LAGUERRE-GAUSSIAN BEAM 


The factor Lp is the associated Laguerre polynomial while 
floo = with P standing for the power of 

the laser beam, L is the atomic excited state decay rate and 
Is the saturation intensity for the atomic transition. Since the 
expression (1) is given in the cylindrical coordinates we also 
have r = y/x^ + and (j) = arctan(?//a:). 

First, we show how we can achieve different phases being 
experienced by two atoms located at different positions. Con¬ 
sider the example shown in Fig. 1 in which a system of two 
atoms is irradiated by a LG beam propagating along the z- 
direction. The shaded ringlike region shows schematically the 
spatial intensity distribution of the beam in the case of p = 0. 



FIG. 1. (Color online) A Laguerre-Gaussian laser beam with a 
doughnut-like intensity profile interacting with two atoms. The beam 
travels along the 2 -axis but is displaced by a distance d along the x- 
axis. The atoms distant ri 2 from each other are located at positions 
(0, ri 2 l 2 , 0) and (0, —ri 2 l 2 , 0), respectively. 


We proceed now to show how we can create a scheme 
where two atoms (emitters) although irradiated by the same 
laser beam, may experience different phases and/or ampli¬ 
tudes of the beam. As mentioned above this can be achieved 
if the driving field has a structured phase and amplitude, i.e. 
a phase and amplitude which have a complex spatial depen¬ 
dence as in the case of an optical vortex beam like an LG or 
a Bessel beam [23]. We choose to work with an LG beam 
whose interaction with a two-level atom result in a Rabi fre¬ 
quency given by [24] 

flop j p\ I rV2 \' ' 
y/l +im+PV- \wiz)J 
X {2r^/w^{z)) (1) 

where we have assumed that the beam propagates in the z- 
direction, w{z) = Wq^/I + z'^jz'j^, wq is the beam waist, 
and r is the position of the atom in the xy plane. The num¬ 
bers I and p are the mode indices with I associated with the 
quantized angular momentum IH carried by the beam photon 
along the propagation axis and p is the radial index associ¬ 
ated with the number of intensity rings in the transverse plane. 


Assume that the interatomic axis has been displaced at a 
distance d along the a;-axis with respect to the center of an LG 
beam and the atoms are located at positions (0, ri 2 / 2 ,0) and 
(0, —ri2/2, 0). The Rabi frequencies fly (j = 1,2) experi¬ 
enced by the atoms can be easily determined from the general 
expression (1), and are given by 
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where j = 1,2. Clearly, there is a phase difference between 
the Rabi frequencies which is due to the azimuthal phase fac¬ 
tors, and is given by 


A(f> = 21 arctan (ri2/2d). (3) 


The phase difference depends on the distance ri 2 between 
the atoms, the lateral displacement d, and the helicity 1. By 
changing one of these parameters we can change the phase 
difference at will. It is easily seen that if there is no displace¬ 
ment along the cc-axis, d = 0, and then the phase difference 
is equal to Acj) = Itt. Of course in a beam with a structured 
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amplitude any change in these parameters will also result in a 
change in the magnitude of the Rabi frequencies. If we keep 
the beam propagation axis symmetrical with respect to the two 
atoms the field amplitudes experienced by the atoms will be 
the same. 

The doughnut-like transverse intensity profile of the LG 
beam has a maximum at radial distances r = wo^/\l\/2, 
where wq is the beam waist at z = 0. The size of the beam 
waist can range from a few hundreds of microns down to half 
the wavelength of the laser beam. To create a phase difference 
equal to tt the angle subtended by the two atoms has to be 
equal to 7r/|^|. In this case the distance between the atoms is 
given by ri 2 = woy^2\l\ sin(7r/2|Z|). For very large values of 
the helicity I, the distance between the atoms is approximated 
by ri 2 = wqtt/^ y2\l\. As we have shown in Ref. [21], for 
very short distances between the atoms their mutual interac¬ 
tion is very crucial for mode switching and routing. These can 
occur when the distance between the atoms is of the order of 
the wavelength of the laser radiation since in this case the in¬ 
teratomic dipole-dipole interaction is significant. Thus if, for 
example, we wish the distance between the atoms to be equal 
to A then the beam waist has to be wq = Xy/2\l\/TT. For a 
value of I = 100, which is experimentally achievable [25], we 
get a beam waist of around 4.5A. Higher values of I can lead 
to larger beam waists. In Ref. [25], LG beams were produced 
with I — 300. The use of such a beam with a waist equal to 
7.8A will ensure that two atoms separated by a distance equal 
to A experience Rabi frequencies with a phase difference equal 
to TT. 

We now demonstrate how one can achieve different Rabi 
frequencies at the positions of the atoms. Since the mutual 
distance between the atoms has to be of the order of the wave¬ 
length for the system to work as a strongly directional antenna 
we must choose proper light fields with amplitude gradients of 
this order of magnitude. The first such case would be the use 
of a standing wave, which could be created along the y-axis. 
In this case we can arrange the configuration in such a way 
that one atom is at the node of the standing wave, thus its Rabi 
frequency is 0, while the other one is at the anti-node of the 
standing wave so its Rabi frequency is the maximum possible. 
This configuration may have several restrictions in its opera¬ 
tion since the distance between the node and the neighboring 
anti-node is always fixed and equal to A/4. This fixed distance 
imposes a restriction on the distance between the atoms. 

The structured beams could give us new opportunities since 
they could have spatial gradients of intensity which, with the 
proper choice of parameters, can be very sharp resulting in the 
desired different Rabi frequencies. The first such case is pro¬ 
vided by an LG beam with a radial number p different from 
0. As noted by Flick et al. [26], this number has been of¬ 
ten overlooked in the literature with the attention of the re¬ 
searchers concentrated on the index I which is related to the 
quantized orbital angular momentum carried by the photon. 
As we know a beam with a radial number p has p -I- 1 max¬ 
ima (rings) in the transverse intensity profile. The distance 
between these rings depends on indices I and p and the choice 
of beam waist. The distance can be analytically calculated for 
generic I and p = 1,2. As our numerical analysis has shown 


the intensity rings come closer as the value of the radial num¬ 
ber p increases while simultaneously we keep I at the lowest 
possible value i.e. I = 1. We chose here the case, where 
p = 20 and 1 = 1. In this case the intensity pattern on the 
transverse plane has 21 rings. In Fig. 2 we show the inten¬ 
sity versus the radial distance from the beam axis. The axis is 
scaled in units of the beam waist wq. It is easily seen that the 
intensity maxima are very close to each other. The first and 
the second maxima are at a distance equal to 0.27u>o. If we 
choose a value for the beam waist equal to 3.7A then the dis¬ 
tance between these maxima becomes equal to the wavelength 
of light A. By placing the two atoms at these points we ensure 
that the Rabi frequency of the first atom can be three times 
larger than that of the second atom. By applying an LG beam 
with a higher values of the radial index the intensity maxima 
can come even closer so we can ensure different Rabi frequen¬ 
cies for atoms at a distance of the order of the wavelength for 
larger values of the beam waist. 
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FIG. 2. (Color online) Intensity pattern (in arbitrary units), up to a 
radial distance equal to the beam waist, in the case of a Laguerre- 
Gaussian laser beam with 1 = 1 and p = 20. The beam is assumed 
to propagate along the 2 axis. In set: How the two atoms have to be 
irradiated by the beam in order to experience different light intensi¬ 
ties and thus to ensure different Rabi frequencies. 

The second case is obtained when we superpose two similar 
LG laser beams, propagating along the z-direction with oppo¬ 
site helicities, I and —1. In this case we get a beam, which 
has a petal-like transverse intensity pattern with a number of 
21 petals. The intensity pattern of this configuration, known in 
the literature as an optical Ferris wheel [26, 27] is illustrated 
in Fig. 3. The inset in Fig. 3 shows the intensity pattern on 
the first quadrant of the xy plane with such a configuration for 
I = 15. It is shown that the intensity goes through 0 as we 
move azimuthally from one point of maximum intensity to 
the next point of maximum intensity. The angle, which sep¬ 
arates these two points is given, in a generic Ferris scheme, 
as A(j> = 180°/|(|. To get two Rabi frequencies with a ratio, 
for example, equal to 10, i.e. G 1 /H 2 = 10, we can place one 
atom at a point of maximum intensity and the other atom at 
a point at an azimuthal angle 6 = 84.3°/|Z|. The length of 
the arc subtended by this angle, which is also the interatomic 
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distance, is ri 2 = (uio/2)\/|Z|/2 sin(84.3°/2|^|). If we use as 
I = 15, then this angle is about 9 = 5.6° and the length of the 
corresponding arc is ri 2 = O.OGTtUQ. This length determines 
the distance between the two atoms. If we wish the distance to 
be equal to half of the wavelength we have to choose a beam 
waist equal to 7.46A. By using LG beams of a higher helicity 
we can ensure larger values of the beam waist. For example, 
for I = 100 the beam waist has to be equal to 19.23A. 



FIG. 3. (Color online) The two atoms located on the xy plane are 
separated from each other by ri 2 and are irradiated with an optical 
Ferris wheel beam propagating along the 2 -direction with helicity 
equal to I = 15. Atoms are at points where they experience different 
intensities of the field. Inset: Intensity pattern at the first quadrant of 
the xy plane. 

The Ferris wheel configuration has yet another interesting 
property. As has been shown [26, 27], if we change the fre¬ 
quency of one of the two beams comprising the Ferris wheel, 
then the intensity pattern rotates in space at an angular fre¬ 
quency ujp = ioji — a;2)/2|(|. When the pattern rotates the 
atoms will be periodically in regions of different intensities, 
thus their Rabi frequencies will vary in time, and this will re¬ 
sult in a case where the intensity pattern emitted by the two 
atoms will rotate in space. 

We should point out that the regime of very small beam 
waists, which we used in our numerical analysis, is associ¬ 
ated with tightly focused light beams carrying orbital angular 
momentum. The use of such beams may introduce effects as¬ 
sociated with the non-paraxial regime [28], which we do not 
take into account in this paper. 


III. GENERAL EORMALISM 


We consider a system composed of two identical emitters 
located at fixed positions, at a distance ri 2 from each other, 
irradiated by a coherent laser beam and interacting with the 
electromagnetic field. Each emitter may be simply a single 
atom, for example, or a quantum dot, or a larger object such 
as a dimer. If atoms are involved, this system can be real¬ 
ized in practice by optical methods, i.e. by considering that 
the atoms are placed at neighboring optical lattice sites or in 


microscopic optical traps [29]. Alternatively, we can consider 
two ions trapped by electromagnetic fields [30]. The latter 
method has played an important role in the demonstration 
of interference effects in the light emitted by two ions [31]. 
Small fluctuations of the atom positions around the trap min¬ 
ima may obscure the effects of mode switching and light rout¬ 
ing but we do not consider them in this work. If we wish to 
avoid the effects of spatial fluctuations, we can consider our 
atoms as generic two-level emitters (for example specially tai¬ 
lored nanoparticles) embedded in a material [20] . We focus on 
a single dipole transition between two non-degenerate energy 
levels of each emitter, the excited je^) and ground \gi) levels, 
and refer to it simply as a two-level system. Thus, the emit¬ 
ters, according to the above description, can be modelled as 
dipoles. Our main concern here is on how the two emitters 
located close to each other and collectively interacting with 
the electromagnetic field can operate as a highly directional 
nanoantenna. We focus on the directional properties of the 
radiation field emitted by the system. 


A. Intensity of the radiation field 


The quantity of central interest is the intensity of the radi¬ 
ation field emitted by two atoms and detected in the far field 
zone of the system. The intensity can be written in terms of 
the positive and negative frequency parts of the electric field as 

^) = t) ■ {R: t)), (4) 

ztt/co 

where fcp = ujq jc, R is a vector pointing in the direction of 
the detection of the field and R is the distance between the 
radiating system and the detector. Here, we have introduced 
the factor (i?^/27rfco) so that /(i?, t)dndt is the probability of 
detecting a photon at time t inside the solid angle element dH. 
around the direction R in the time interval dt. 

The electric field radiated by the atoms can be expressed in 
terms of the atomic dipole moments. The negative and pos¬ 
itive frequency parts of the field in the far field zone of the 
radiating atoms can be written as the sum of dipole fields. 


E'^-'>{R,t) 





2=1 


[Ri X {Ri X jli)] 

r3 e 

[Ri X (Ri X fii)] 


(5) 


where fli is the transition dipole moment, Sf and S~ are the 
usual raising and lowering operators of atom i, Ri is the posi¬ 
tion vector of atom i, and R is the unit vector in the direction 
of observation [R = R/R). 

When Eq. (5) is substituted into Eq. (4), the intensity of 
the radiation field measured in the direction R at time t be¬ 
comes [32, 33] 


2 

I{R,t) = u{R) ^ (5+(f)5-(f))e“-- , 


*j'=i 


( 6 ) 
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where u{R) = (SF /Sir) sin^ r?, with 1 } the angle between the 
direction of observation R and the direction of the atomic 
dipole moments, and F— the spontaneous emission (damp¬ 
ing) rate of the atomic transition. It is seen that the intensity 
of the radiation field is determined by the correlation functions 
of the atomic dipole operators. The intensity can be written in 
terms of four contributions involving the correlation functions 
and geometrical factors 

liRR) = u{R){{Stit)S^it)) + {Stit)S^it)) 

+ {t)S 2 (t)) + {S 2 {t)S^ (f))] cos(fcri 2 cos 9) 

+ i[{S^{t)S2it))-{S:^{t)S^{t))]sm{kri2COs6)} , (7) 

where 9 is the angle between fi 2 and the direction of obser¬ 
vation R. The variation of the intensity with the observation 
angle 9 is called the radiation pattern. Certain general fea¬ 
tures of the radiation pattern follow from Eq. (7). The first 
term expresses the intensity of the emitted radiation created 
by spontaneous emission from atom 1. The second term ex¬ 
presses the intensity of the radiation spontaneously emitted 
by atom 2. These two terms are always positive and indepen¬ 
dent of 9. The contribution of these two terms obviously leads 
to a spherical shape of the radiation pattern. The third and 
fourth terms result from the interference between the radia¬ 
tion fields emitted by different atoms. If nonzero, these terms 
can lead to a non-spherical shape of the radiation pattern and 
the emitted radiation can exhibit a strong enhancement or re¬ 
duction in a direction 9 at which cos(A:ri 2 cos 9) = ±1 and/or 
sin(A:ri 2 cos 0) = ±1. 

Let us look at some features of the cos(fcri 2 cos 9) and 
sin(A:ri 2 cos 0) factors in I{R,t) which define directions of 
maximum and minimum emission. First, since the cosine and 
sine functions are shifted in phase by 7r/2, we see that the di¬ 
rections in which these two terms can enhance or reduce the 
intensity do not overlap. In particular, for atoms separated 
by a distance ri 2 = A/2, the factor cos(A:ri 2 cos 0) = 1 for 
0 equal to tt/2 and 3tt/2, while sin(A:ri 2 cos 0) = 1 for 0 
equal to 7r/3 and Stt/S. Moreover, the number of directions 
in which the intensity can be enhanced or reduced increases 
with an increasing ri 2 . For instance, when ri 2 = A, the fac¬ 
tor cos(fcri 2 COS0) = 1 for 0 equal to 0,7r/2, tt, and 37r/2, 
whereas sin(fcri 2 cos0) = 1 for 0 equal to 0.427r and I.SStt. 
Thus, if the goal is for the system to work as a highly direc¬ 
tional nano-antenna emitting light only in a few directions, 
then the emitters should be kept at distances shorter than the 
resonant wavelength, ri 2 < A. 

There is an another important difference in the directional 
properties of the two factors. It is easy see that if there is a 
direction 0 in which cos(fcri 2 cos0) is maximal (= 1), it is 
also maximal in the opposite direction 0 -f 180°. This means 
that the factor cos(/cri 2 cos 0) has the property of concentrat¬ 
ing the radiation along two opposite axial modes. We refer to 
this feature as an axial concentration of the radiation or a two- 
sided emission. The directional property of sin(fcri 2 cos 0) is 
different. If there is a direction 0 in which sin(fcri 2 cos 0) is 
maximal, it is also maximal in the direction 9' = 360° — 0. 
Thus, the factor sin(fcri 2 cos 0) has the property of concen¬ 
trating the emission on one side of the system in two axial 


modes propagating in directions differing by 20. We refer 
to this feature as a spatial concentration of the radiation or a 
one-sided emission. The differences in the ways in which the 
factors cos(fcri 2 cos 0) and sin(fcri 2 cos 0) affect the radiation 
pattern are illustrated graphically in Sec. IV. 

It is convenient, in particular for physical interpretation, to 
write the intensity in terms of the density matrix elements of 
the density operator of the two-atom system represented in the 
basis of the superposition (collective) states [34-36] 

\g) = \ 9 i)\ 92 ), |e) = |ei)|e 2 ), 

k) = ^(|ei)|02) + |ffi)|e2)), 

|a) = ;^(|ei)|02)-|5i)|e2)), (8) 

where |s) and |a) are, respectively, the symmetric and anti¬ 
symmetric combinations of the product of bare atomic states. 
Here, je^) and \gi) represent the excited and ground states of 
atom i. It is easily shown that in terms of the density matrix 
elements the correlation functions appearing in Eq. (7) are 

{S^{t)S^{t))-\-{S}{t)S 2 it)) = Pssit)-\-paa{t)-\-2pee{t), 

(S'/'(f)5'2"(f)) + (S'^(f)5'r (t)) = p,sit) - Paa(t), 

I [{St(t)S 2 {t)) - {St(t)S^(t))] = 21m[pas{t)]. (9) 

We see from Eq. (9) that the interference term proportional 
to cos(fcri 2 COS0) will contribute to the intensity only when 
Pssit) Paa(,t), i-C- when the symmetric and antisymmet¬ 
ric states are unequally populated. Consequently, a reduc¬ 
tion in the population of one of the two states relative to the 
population of the other state will be accompanied by an ax¬ 
ial two-sided emission. On the other hand, the interference 
term proportional to sin(fcri 2 cos0) will contribute to the in¬ 
tensity only when Im[pas(f)] 7 ^ 0. Thus, a nonzero coherence 
between the symmetric and the antisymmetric states will be 
accompanied by a one-sided emission. 

B. Master equation 

Suppose that the emitters are located at positions, fi = 
(0, —ri 2 / 2 , 0) and r 2 = ( 0 ,ri 2 / 2 , 0 ) along the y axis and 
are illuminated with a monochromatic laser beam of angular 
frequency ujl, a propagation wave vector Pl, and of a variable 
intensity profile with spatially varying amplitude El{x, y, z) 
and phase 4>l{x, y, z), as shown in Fig. 4. In areal experiment 
this might be an LG laser beam. The laser excites transitions 
between the two energy states je^) and \gi) {i = 1,2), sepa¬ 
rated by a frequency ujq. The population of the excited states 
of the emitters decays to the ground states at a rate F. 

The emitters radiate by the process of spontaneous emis¬ 
sion and the intensity of the emitted field is detected at a point 
A distance R in the yz plane. To explore the role of the struc¬ 
tured beam profile in the controlled mode switching and direc¬ 
tional emission, we choose the direction of propagation of the 
laser beam perpendicular to the interatomic axis, _L fi 2 . 
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FIG. 4. (Color online) Schematic diagram of the system composed of 
two emitters represented hy transition dipole moments /2i and pL 2 and 
irradiated with a laser beam propagating in the direction perpendic¬ 
ular to the interatomic axis, Isl -L fi 2 . The laser has a cylindrically 
symmetric (doughnut shaped) intensity profile. With this variable 
intensity profile, the two emitters can be excited with arbitrary am¬ 
plitudes and phases. 


i.e. the laser beam propagates along the z-direction and has 
a structured amplitude in the xy plane. If the driving field 
is in the form of a plane wave of a constant amplitude then 
the atoms experience the same amplitude and phase of the 
field. However, when the laser field is taken to have a vari¬ 
able intensity profile, then the atoms can experience different 
amplitudes and phases of the field; as a result, the radiative 
properties of the atoms may change. 

Our purpose is to calculate the angular distribution of the 
radiation intensity, the radiation pattern of the field emitted 
by the driven atoms. Thus, according to Eq. (7), we need to 
obtain the atomic correlation functions or the matrix elements 
of the reduced (atomic) density operator, which satisfies the 
master equation [32, 33] 


dp 

dt 


i 

h 


[Hq + Hl, p] + 




( 10 ) 


where 


Hq — HAl {Ass+ Aaa+‘2'Aee) + fi^l2 ~ ^aa) j (H) 


and 


Hl = 


n 


{Hi (Ags + Asg T Aag ^ea) 


2V2 

+ H 2 (Ags + Agg — Aag + Aea) + H.C.} 


( 12 ) 


is the interaction Hamiltonian of the atoms with the driving 
laser field. Here, Anm = W) are projection operator be¬ 
tween the collective states, = ojq — ujl is, the detuning 


of the laser field frequency lol from the atomic transition fre¬ 
quency Wo, and the quantities Hi and H 2 are the Rabi fre¬ 
quencies of the laser field at the position of the atom 1 and 
2, respectively. The quantity H 12 depends on the distance be¬ 
tween the atoms and gives the effect of the atomic interaction, 
the dipole-dipole interaction, on the shift of the energy levels 
of the system 


Hi2 — -T ■( — 


1 - (m ■ ri2) 


cos{kori2) 


1 - 3 (/r • ri 2 )^ 


kori2 

sin(fcori2) cos(fcori2) 
{kori2f {kori2f 


(13) 


where fcp = Wq/c, and /2 and ri2 are unit vectors in the direc¬ 
tion of the atomic dipole moment fl and the interatomic axis 
ri2, respectively. 

The dissipative part of the master equation (10) consists of 
two terms corresponding to the two decay (fluorescent emis¬ 
sion) channels [32, 34]; the symmetric channel |e) —>■ |s) —>■ 
\g) with an enhanced decay rate Tg = T + ri 2 : 

+P {Aee + ^ss) + AggpAgg + AgspAes 

~‘^{AsePAes -i- AgspAgg)} , (14) 

and the anti-symmetric channel |e) —>■ |a) — \g) with a re¬ 
duced decay rate = T — ri 2 : 

~ {(^ee + Aaa) P 

+p ee + A 

aa ) ^aeP-^ag -^gaP-^ ea 

2 (AaepAgcj + AgapAag^y . (15) 


Here, ri 2 gives the effect of the atomic interaction on the 
damping rate of the system 


ri2 = 


1 - (p ■ ^12) 


sin(A:o?’i2) 


1 - 3(^ • ri2)^ 


kori2 

cos(fcori2) sin(A:ori2) 


ikori2y 


{kori2)'' 


(16) 


Since the Rabi frequencies Hi = |Hi| ex-p{i(j)i) and H 2 = 
IH 2 I ex-p{i(j) 2 ) can be different, i.e. can have different magni¬ 
tudes and phases, we write the interaction Hamiltonian (12) in 
the form 


Efi — ik [Hc[(Aes+Asg) + H^(Aea+A(jg) — H.C.] , (17) 


where 


H(j — (Hd sin (l)fi zHq cos j 

H/3 = (Hosin(;)d-*i2dCOS(;)d)/V2, (18) 
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in which 


^0—2 (l^il + 1 ^ 2 !)) 2 (l^il ~ 1 ^ 2 !), (19) 


and (f>d = {(f>i — 4 > 2)/2 is the phase difference between the 
Rabi frequencies of the atoms. Clearly, the transitions of the 
symmetric channel are driven at the Rabi frequency while 
the transitions of the antisymmetric channel are driven at the 
Rabi frequency In general, both channels can be simulta¬ 
neously driven by the laser. 


C. Equations of motion for the density matrix elements 


In the basis of the collective states the master equation leads 
to a set of 15 coupled differential equations for the density 
matrix elements. Among them we can distinguish eight equa¬ 
tions for the density matrix elements determining transitions 
between the symmetric states 


Pss - 

Pee — 

Pag = 

Pse — 

Peg = 
Pgs = 


Pea — 


Pee') ^^a)Pga Pae) “f C.C.] , 

2Tp ee “t“ [^a/^se - ^pp ae C-C-] ■> 

^ 2 ^*^ ^^12^ Psg ~^^sPes H” ^aiPgg Pss') 

^aPeg ^0P sa 5 


- (2rs + Ta) + i^l2 


Pse ^aiPss pee) 


Pge ^pP sa: 

^Peg ^a:(Psg Pes ) - ^p(p ag + p ea) •) 

^2^*^ Pgs sPse ^OL^Pgg Pss) 

Pge ^pp as 5 

1 


(2rs + Ta) — iTti2 


Pes 'T ^aiPss Pee) 


^aPeg ^0P as 1 

Pge — ^Pge “t“ ^a^Pgs Pae) ^piPga Pae)^ (20) 


and seven equations determining transitions between antisym¬ 
metric states 


Paa 

Pae 


^a{^Paa Pee) ~T “1“ Pae) ~T C.C.] 


, (rs+2ro) — iTti2 


“h Tl^Pg^ “h Tl^Pas-) 


Pae ^pi^Paa Pee) 


Pag — (^ 2 ^^ ^^ 12 ^ Pag ^ 2 ^aPea ~T Ttpi^Pgg Paa) 

-\- Tl^p^g ^aPas 5 

Pas — (r 2 zrii 2 )pas ^cx.Pae ~T Tl^Pag 

+ TlpPgs + 

es j 


Pea — 


-(rs-i-2ro) + iVii2 


Pea ^pi^Paa Pee) 


Pga — 


Paa — 


-f ^^Peg -f ^aPaa^ 

^2^^ ^1^12^ Pga aPae “t” ^piPgg Paa) 

^pPge ^aPaai 

(r “t“ 2i^i2)paa ^oiPea “f ^ctPga 

+ ^*pPag + ^PPae- (21) 


The remaining equation of motion for pgg is found from the 
closure relation pgg = 1 — psa — Pee — Paa- We see from 
the equations that the transitions between the symmetric states 
are driven at the Rabi frequency and are coupled to the 
antisymmetric states with a strength proportional to 17^. On 
the other hand, transitions between the antisymmetric modes 
are driven at the Rabi frequency Q.p and are also coupled by 
to the transitions between the symmetric states. Thus, in 
the case of = 0, the dynamics of the symmetric and anti¬ 
symmetric states are independent of each other. 

For numerical analysis, it is convenient to write the set of 
differential equations in matrix form. When pgg is eliminated 
from the equations, we arrive at an inhomogeneous equation 


—Y = -MY + P, (22) 

at 

where F is a column vector composed of the density matrix 
elements 


F — COl(pss, Pee j Psg: Pse ? Peg ? Pgs j Pes i Pge i 

Paai Paei Pag-} Pasi Pea-) Pga-) Psa ), ( 23 ) 

P is a column vector with nonzero elements 

P 3 = Pe = nl, Pn = P 14 = ( 24 ) 

and M is the 15 x 15 matrix of the complex coefficients. It is 
convenient to express the matrix M in block form as 

M=(*®). (25) 

in which the block S' is an 8 x 8 matrix of the coefficients 
of the eight equations involving the density matrix elements 
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(20), block A is a 7 X 7 matrix of the coefficients of the seven 
equations involving the density matrix elements (21). 

The off-diagonal blocks B and D are 9x6 and 6x9 ma¬ 
trices whose nonzero elements are 

B 21 = .B 37 = i?53 = i?55 = i?74 = rifj, 

B 25 = i?47 = ^64 = Bs 2 = Bsq = (26) 

and 

Di2 = —Ta, 

D28 = = Dqs = -D74 = —n/ 3 , 

D31 = D32 = D32 = ^/ 3 , 

D22 = —D33 = £>55 = —£>61 = —Dq2 = d^3 = — 

(27) 

Clearly, the nonzero elements of the off-diagonal blocks are 
only those involving Vlp and its complex conjugate. Note that 
Up is different from 0 only when the atoms experience differ¬ 
ent amplitudes and/or phases of the driving laser field. 

IV. RADIATION PATTERN 

We now proceed to give illustrative examples of the direc¬ 
tionality of the emission by a system of two identical two-level 
emitters driven by an LG beam. For this purpose, we numer¬ 
ically solve the set of the optical Bloch equations, Eq. (22) 
to obtain the steady-state {t —)■ 00 ) values of the density ma¬ 
trix elements. The solutions are then applied to determine the 
steady-state radiation intensity. Written in spherical coordi¬ 
nates the radiation pattern of a two-atom system, Eq. (7), is 
independent of the azimuthal angle (j). Consequently, the ra¬ 
diation pattern is displayed graphically in a polar form for a 
variety of different sizes of the system, i.e. different distances 
between the atoms. It can be regarded as the radiation pattern 
viewed from the top of the two-atom antenna. The following 
interatomic distances are chosen for specihc consideration, 
?'i 2 = A/4, A/2,3A/4 and A. The choice of short distances 
(ri 2 < A) has been dictated by the fact that the directional 
character of the emitted field results from the interference be¬ 
tween the electric fields emitted by the different atoms which, 
on the other hand, is pronounced for small ri 2 . Moreover, at 
distances ri 2 < A the system emits light only in a few discrete 
directions, which justifies regarding the systems as a highly 
directional antenna. Radiation patterns are illustrated for two 
particularly interesting configurations of the driving field at 
the position of the atoms, i.e. when the driving held creates 
a phase difference only 7 ^ 0,12^ = 0), and an amplitude 
difference only ( 12 ^ ^ 0 , 4 >d = 0 ). 

Considering the phase difference hrst, we plot in Eig. 5 the 
angular distribution of the radiation intensity for the atomic 
separation ri 2 = A/4 and phase differences, = 0 and 
(j)d = ±7r/4. The case of (pd = 0 corresponds to the dipole 
moments of the atoms driven with the same phase. The an¬ 
gular distribution of the emitted radiation is asymmetric with 
an enhanced emission in the direction perpendicular to the in¬ 
teratomic axis. Thus, for pd = 0, the system tends to radiate 


90 



FIG. 5. (Color online) Polar diagram of the radiation intensity as 
a function of the observation direction 6 for the atomic separation 
ri 2 = A/4 with Al = 0, fii = G 2 = 0.2r and different val¬ 
ues of the phase difference pd- Pd = 0 (dashed-dotted black line), 
= 7r/4 (solid red line), and pd = —7r/4 (dashed blue line). The 
positions of the atoms with the transition dipole moments polarized 
perpendicular to the plane of the paper {p, _L fi 2 ) are shown by filled 
(brown) circles. 


in the direction perpendicular to the interatomic axis. In other 
words, the system behaves as a two-sided antenna radiating 
along the direction perpendicular to the interatomic axis. Note 
that the pattern is symmetric around both the horizontal (0, tt) 
and the vertical {tt/2, 3tt/2) axis. The pattern varies with the 
phase difference pd with which the atoms are driven. We see 
that for the phase differences pd = ±7r/4 the system tends to 
radiate along the interatomic axis with a striking asymmetry 
about the vertical axis. Depending on whether pd = 7r/4 or 
pd = — 7r/4, the emitted radiation is spatially concentrated in 
the left half (cos 0 < 0) or in the right half (cos 0 > 0) of 
the pattern. Thus, the direction of the emission reverses when 
pd reverses sign from 7r/4 to —7r/4. Clearly, a phase differ¬ 
ence between the atomic dipole moments dictates the direc¬ 
tion of the emission. Eor pd = ±7r/4 the system behaves as 
an one-sided antenna. Another interesting observation is that 
the radiation patterns for pd = ±7r/4 resemble very much the 
radiation pattern of a Yagi-Uda antenna [37]. Similar pattern 
is also produced by a patch antenna [38]. 

The switching of the behavior of the system from a two- 
sided towards a one-sided antenna can be understood by in¬ 
voking the effect of the angular factors cos(/cri 2 cos 9) and 
sin(fcri 2 cos 6) that determine the radiation pattern, Eq. (7). 
As discussed in Sec. Ill A, the factor cos(A:ri 2 cos0) has the 
property of concentrating the radiation along two opposite ax¬ 
ial modes, whereas the factor sin(A:ri 2 cos 9) has the property 
of concentrating the emission along single axial modes. Ac¬ 
cording to Eqs. (7) and (9), the factor cos(A:ri 2 cos 9) con¬ 
tributes to the radiation pattern only when the symmetric and 
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antisymmetric states of the system are unequally populated, 
Pss — Paa 7^ 0. The factor sin(fcri 2 cos 9) contributes to 
the radiation pattern only when there is a nonzero coherence 
between these states, Im[pos(i)] 7^ 0. It is easy to find 
that for the parameter values in Fig. 5, pss — Paa 7 ^ 0 for 
both (j)d = 0 and (pd — ±7r/4, whereas Im[/3as(f)] = 0 
for (pd = 0, but Im[pQ 5 (t)] ^ 0 for pd = ±7r/4 such 
that Im[pas(f)]|0rf=7r/4 = -Ini[pas(t)]U,i=-,r/4- Thus, the 
switching of the behavior of the system from a two-sided 
towards a one-sided antenna is achieved by the creation of 
nonzero coherence between the symmetric and antisymmetric 
states. 

In Fig. 6 we compare the angular distribution for pd = 0 
with that for pd = 7r/2 for the same atomic separation as 
in Fig. 5. The phase difference pd = 7r/2 corresponds to 
the dipole moments of the atoms driven with opposite phases, 
pi = TT, (^2 = 0. As discussed in the previous section, the case 
of pd = 0 corresponds to the situation where the laser field 
drives only the symmetric modes of the system, while in the 
case of pd = tt/2 the laser effectively drives only the antisym¬ 
metric modes. We see that the phase difference pd = tt/2 has 
the effect of forcing the atoms to radiate along the interatomic 
axis. Unlike the radiation patterns for the phase differences 
pd = ±7r/4, this radiation pattern exhibits highly pronounced 
lobes along the interatomic axis with zeros (nodes) in the 
emission occurring in the perpendicular directions, 6* = 7r/2 
and 6 — 37r/2. The opening angle of the lobes is about 120°, 
which is much smaller than that for pd = ±7r/4. Behavior 
of this kind can be interpreted as being a consequence of the 
trapping of photons along the interatomic axis: a photon emit¬ 
ted by atom 1 is absorbed by atom 2 and then when emitted 
by atom 2 is re-absorbed by atom 1. 

90 



FIG. 6. (Color online) Polar diagram of the radiation intensity as 
a function of the observation direction 0 for the atomic separation 
ri2 = A/4 with = 0, = U2 = 0.2r and different phase 

differences pd'- pd = Q (dashed-dotted black line), pd = 7r/2 (solid 
red line). 

We would like to stress here that although the emission 


for pd = 1^/2 is more focused along the interatomic axis than 
that for (/)d = 7r/4, it is symmetrically distributed around the 
vertical axis. For this reason, the directionality of the emis¬ 
sion for pd = 7r/4, which is strongly asymmetric around the 
vertical axis, can still be regarded as more pronounced than 
that for pd = 7r/2. 

90 



FIG. 7. (Color online) Polar diagram of the radiation intensity as 
a function of the observation direction 0 for the atomic separation 
ri 2 = A/2 with = 0, Gi = G 2 = 0.2r and for different phase 
differences pd’- pd = (solid black line), pd = 7r/4 (dashed-dotted 
red line), and pd = it 12 (dashed blue line). 

The directionality of the emission by the two-atom sys¬ 
tem can be improved by increasing the separation between 
the atoms. This is illustrated in Fig. 7 which shows the radia¬ 
tion pattern for the atomic separation ri 2 = A/2 and various 
values of pd- An improvement of the directionality shows up 
clearly in the presence of zeros in the emission for several 
phase differences pd- The radiation patterns are composed of 
two well distinguished lobes whose directions and symmetry 
depend on pd- A change of pd from 0 to 7r/2 rotates the direc¬ 
tion of the lobes by 7r/2; from the perpendicular to the parallel 
to the interatomic axis. For phases 0 and tt/2 the lobes are the 
same size but for any phase between these two values they are 
different sizes. Moreover, the opening angles of the lobes and 
their magnitudes also vary with pd- 

At a larger distance between the atoms a variation of the 
phase difference pd not only may lead to a change in the di¬ 
rection of the lobes but also may reduce the number of lobes. 
This is illustrated in Fig. 8, which shows the radiation pattern 
for ri 2 = 3A/4 and two values of pd- It is shown that the 
number of lobes depends on pd and is reduced from four to 
three when pd is varied from 0 to 7r/4. Moreover, we see that 
upon changing the phase pd from 0 to 7r/4, the behavior of 
the system changes from that of a two-sided antenna to that of 
a one-sided antenna. 

However, the reduction in the number of lobes seen for 
ri 2 = 3A/4 might not be seen for other distances between 
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FIG. 8. (Color online) Polar diagram of the radiation intensity as 
a function of the observation direction 6 for the atomic separation 
ri 2 = 3A/4 with A/. = 0, = fl 2 = 0.2r and different phase 

differences cfyd'. 4>d ~ 0 (dashed black line), (pd = rr/d (solid red 
line). 

the emitters. For example, at distance ri 2 = A, the varia¬ 
tion of the phase (j)d from 0 to 7r/2 changes the direction of 
the lobes and their shape, but the number remains the same. 
This is illustrated in Fig. 9 for various values of (pd from 0 
to 7r/2. The figure shows that for pd equal to 0 and 7r/2 the 
radiation pattern exhibits pronounced lobes in directions sym¬ 
metrically located about the horizontal and vertical axis. The 
pattern becomes asymmetric for values of pd different from 0 
and tt/2. Consequently, for pd = 0, tv/2, the system behaves 
as a strongly directional two-sided antenna with the emission 
into narrow lobes. A change in phase from those of 0 or 7r/2 
causes the system to behave as an one-sided antenna. 

The results presented in Figs. 5—9 show clearly that in 
all cases considered it turns out that the system behaves as 
a two-sided directional antenna when the phase differences 
(pd equals an even multiple of 7r/4. When (j>d equals an odd 
multiple of 7r/4, the system behaves as an one-sided antenna. 
It is particularly simple to interpret this behavior by referring 
to the expression for the radiation intensity, Eq. (7). The in¬ 
tensity contains two interference terms, one proportional to 
cos(fcri 2 cos 9) and the other proportional to sin(fcri 2 cos 9). 
Radiation at an angle 9 from atom 2 has a distance ri 2 cos 9 
farther to travel to the observation point than that from atom 
1. Its electric field is therefore retarded by an extra amount 
kr 12 cos 9. The fields of the atoms interfere constructively if 
the phase (/)2 at which the atom 2 is driven relative to phase pi 
at which atom 1 is driven, pi — p 2 = ‘2-pd, can compensate 
for the additional retardation. Thus, constructive interference 
occurs at angles 9 such that 

kri 2 cos 9 = 2pd. (28) 

For the phase difference pd = 0 or 7r/2, the interference 


90 



90 



FIG. 9. (Color online) Polar diagram of the radiation intensity as 
a function of the observation direction 9 for the atomic separation 
ri 2 = A with Al = 0, fli = 0.2 = 0.2r, and different phase 
differences pd- Top: The dashed black line represents pd = 0, the 
solid (red) line pd = tx 12. Bottom: The dashed black line represents 
pd = 7r/8, the solid (red) line pd = ix/4.. 


term cos(fcri 2 cos 0) = ±1 whereas sin(fcri 2 cos 0) = 0. 
In this case, the variation of the radiation intensity with 9 
is solely determined by the cosine term. Thus, a two-sided 
emission takes place simply because cos 0 > 0 gives the same 
cos(fcri 2 COS0) as cos0 < 0. However, for the phase dif¬ 
ference pd = 7r/4, the interference term cos(A:ri 2 cos 0) = 
0 whereas sin(fcri 2 cos 6*) = 1. Hence, the symmetry of 
the emission direction is broken simply because cos 9 > 
0 gives sin(A:ri 2 cos 9) positive whereas cos 9 < 0 gives 
sin(A:ri 2 cos 9) negative. The radiation intensity then has dif¬ 
ferent magnitudes in the cos 0 > 0 and cos 9 < 0 parts of 
the pattern. In physical terms one-sided emission reflects con¬ 
structive and destructive interference effects on the two sides 
of the pattern. 

The above analysis of the radiation pattern has been fo- 
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FIG. 10. (Color online) Polar diagram of the radiation intensity as 
a function of the observation direction 6 for the atomic separation 
ri 2 = A/4 with A_l = 0 and i/i = (/>2 = 0 and for different Rabi 
frequencies experienced by the atoms, fli = O.dF, fl 2 = 0 (dashed 
black line), fli = 0, G 2 = O.dF [solid (red) line]. 

cused on examples where the emitters experience the same 
amplitudes but different phases of the driving field. Inter¬ 
esting behavior can be uncovered in the radiation pattern by 
considering a situation where the atoms experience different 
amplitudes of the driving field, Q,d 7 ^ 0. We illustrate this 
with an example in which the driving of the system is config¬ 
ured such that only one of the two atoms is driven by the laser 
field. Referring to the collective modes of the system discus¬ 
sion in Sec. IIIC, this corresponds to the situation where the 
symmetric and antisymmetric modes are driven by the same 
Rabi frequency, In Fig. 10 we show the radiation 

pattern for ri 2 = A/4, = 0. The radiation pattern is com¬ 

posed of a single lobe showing that the radiation tends to be 
mainly on one side of the system. Thus, the un-driven atom 
steers the system to radiate towards one side of the pattern, 
i.e. it shows a tendency to behave as an one-sided antenna. 
The direction of the lobe reverses when 17^ reverses its sign. 
Again, this one-sided emission behavior can be understood 
by referring to the effect of the angular factor sin(A:ri 2 cos 9) 
appearing in the expression for the radiation pattern. The 


contribution of this factor is determined by Im[/ 3 as(f)]. It 
is easily shown that Im[pas(f)] 7 ^ 0 for 17^ 7 ^ 0 such that 
Im[pa,,(7)]|Q^>o = -Im[pa^(7)]|Q^<o. 


V. SUMMARY 

We have studied the directional properties of the radiation 
field emitted by a system of two identical two-level emitters 
located at short distances and driven by a laser beam with 
structured phase and amplitude. We have shown that the 
system can operate as a nanoantenna for controlled directional 
emission. We have calculated the radiation intensity of the 
field emitted by the system and have shown that a constant 
phase or amplitude difference at the positions of the emitters 
plays an essential role in the directivity of the emission. 
Polar diagrams have been presented showing the radiation 
patterns under various conditions of excitation and for various 
separations between the emitters. We have demonstrated 
that depending on the phase or amplitude difference at the 
positions of the emitters, the system can operate as a two- or 
as a one-sided directional nanoantenna. A two-sided direc¬ 
tional emission is achieved with a symmetric driving where 
the emitters experience the same amplitude and a constant 
phase difference of the driving field. We have found a general 
directional property namely that when the phase difference at 
the positions of the emitters equals an even multiple of 7 r/ 4 , 
the system behaves as a two-sided antenna. On the contrary, 
when the phase difference equals an odd multiple of 7 r/ 4 , 
the system behaves as an one-sided antenna. We have also 
considered the case where the emitters experience the same 
phase but different amplitudes of the driving field and have 
found that the effect of different amplitudes is to cause the 
system to behave as an one-sided antenna radiating in one 
direction centered along the interatomic axis. 
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